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Prepotential approach to solvable rational potentials 
and exceptional orthogonal polynomials 

Choon-Lin Ho 

Department of Physics, Tamkang University, Tamsui 25137, Taiwan, R.O.C. 

We show how all the quantal systems related to the exceptional Laguerre and Jacobi 
polynomials can be constructed in a direct and systematic way, without the need of shape 
invariance and Darboux-Crum transformation. Furthermore, the prepotential need not be 
assumed a priori. The prepotential, the deforming function, the potential, the eigenfunctions 
' and eigenvalues are all derived within the same framework. The exceptional polynomials are 

' expressible as a bilinear combination of a deformation function and its derivative. 

§1. Introduction 

\ In the last three years or so one has witnessed some interesting developments 

' in the area of exactly solvable models in quantum mechanics: the number of ex- 

actly solvable shape-invariant models has been greatly increased owing to the dis- 
I covery of new types of orthogonal polynomials, called the exceptional X£ polynomi- 

g^jgflJ-lTSli U]2ii]jg iiiQ classical orthogonal polynomials, these new polynomials have 
1^ . the remarkable properties that they still form complete sets with respect to some 

! positive-definite measure, although they start with degree i polynomials instead of 

a constant. 

Two families of such polynomials, namely, the Laguerre- and Jacobi-type Xi 
•r^ I polynomials, corresponding to £ = 1, were first proposed by Gomez-Ullate et al. 

' in Ref. |T|) , within the Sturm-Lioville theory, as solutions of second-order eigenvalue 

equations with rational coefficients. The results in Ref. JJ were reformulated in 
l/^ I the framework of quantum mechanics and shape-invariant potentials by Quesne et 

■ al.'2J''El' These quantal systems turn out to be rationally extended systems of the 

. traditional ones which are related to the classical orthogonal polynomials. The most 

general Xi exceptional polynomials, valid for all integral £ = 1,2,..., were discovered 
by Odake and SasakP (the case of ^ = 2 was also discussed in Ref. El)- Later, in 
Ref. [5| equivalent but much simpler looking forms of the Laguerre- and Jacobi-type 
Xi polynomials were presented. Such forms facilitate an in-depth study of some 
important properties of the X^ polynomials, such as the actions of the forward and 
. backward shift operators on the X£ polynomials, Gram-Schmidt orthonormalization 

for the algebraic construction of the X^ polynomials, Rodrigues formulas, and the 
generating functions of these new polynomials. Structure of the zeros of the exception 
polynomials was studied in Ref. |6|). 

More recently, these exceptional polynomials have been studied in many ways. 
For instance, possible applications of these new polynomials were considered in 
Ref. Uj) for position-dependent mass systems, and in Ref.iSj) for the Dirac and Fokker- 
Planck equations. The new polynomials were also considered as solutions associated 
with some conditionally exactly solvable potentials.'^ These polynomials were re- 
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cently constructed by means of the Darboux-Crum transformation P'^^^'H'' Rational 
extensions of certain shape-invariant potentials related to the exceptional orthogo- 
nal polynomials were generated by means of Darboux-Backlund transformation in 
Ref. tl2j). Generalizations of exceptional orthogonal polynomials to discrete quantum 
mechanical systems were done in Ref . 1131) . Structure of the X^^ Laguerre polynomials 
was considered within the quantum Hamilton- Jacobi formalism in Ref. I14p . Gener- 
alizations of these new orthogonal polynomials to multi-indexed cases were discussed 
in Ref. I16p J17p . Recently radial oscillator systems related to the exceptional Laguerre 
polynomials have been considered based on higher order supersymmetric quantum 
mechanics.'^ 

So far most of the methods employed to generate the exceptional polynomials 
have invoked in one way or another the idea of shape invariance and/or the related 
Darboux-Crum transformation. Furthermore, the so-called superpotentials (which 
we shall call the prepotential) , which determine the potentials, have to be assumed 
a priori (often with good educated guesses). 

The aim of this paper is to demonstrate that it is possible to generate all the 
quantal systems related to the exceptional Laguerre and Jacobi polynomials by a 
simple constructive approach without the need of shape invariance and the Darboux- 
Crum transformation. The prepotential (hence the potential), eigenfunctions and 
eigenvalues are all derived within the same framework. We call this the prepotential 
approach, which is an extension of the approach we employed to construct all the 
well-known one-dimensional exactly solvable quantum potentials in Ref. I19p . 

The plan of this paper is as follows. Sect. 2 presents the ideas of prepotential 
approach to systems which are rational extensions of the traditional systems related 
to the classical orthogonal polynomials. Li Sect. 3 the prepotential approach is em- 
ployed to generate the L2 Laguerre system. Construction of the Jl and J2 Jacobi 
cases is then outlined in Sect. 4. Sect. 5 summarizes the paper. Appendix A col- 
lects some useful results on the classical Laguerre and Jacobi polynomials. The LI 
Laguerre system is then summarized in Appendix B. 



§2. Prepotential approach 

2.1. Main ideas 

We shall adopt the unit system in which h and the mass m of the particle are 
such that h = 2m = 1. Consider a wave function 0(x) which is written in terms of 
a function W{x) as 

0(2;) = e^W. (2-1) 

Operating on (p^ by the operator —dP/dx^ results in a Schrodinger equation T-Lcf) = 0, 
where 

« = -|, + V-, (2.2) 

V = W'^ + W. (2-3) 
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The dot denotes derivative with respect to x. Since determines the potential 

V, it is therefore called the prepotentialP^ For clarity of presentation, we shall often 
leave out the independent variable of a function if no confusion arises. 
In this work we consider the following form of the prepotential 



W{x, rj) = Wo{x) — ln,^(r/) + \np{r]). 



(2-4) 



Here Wo(a;) is the zero-th order prepotential, and ??(x) is a function of x which 
we shall choose to be one of the sinusoidal coordinates, i.e., coordinates such that 
rj{x)'^ is at most quadratic in ry, since most exactly solvable one-dimensional quantal 
systems involve such coordinates. The choice of 7/(x) and the final form of the 
potential dictate the domain of the variable x. ^{rj) and p{rj) are functions of rj to be 
determined later. We shall assume ^(r/) to be a polynomial in rj. The function p{r]) 
consists of the eigen-polynomial, but itself need not be a polynomial (see Sect. 3). 
With the prepotential (j2-4p . the wave function is 



aWoix) 



■p{ri), 



(2-5) 



and the potential V = + W takes the form 



V = W^ + Wo + 



7?V+ ( 2Wo?7 + ^?-2?7^^ ]p' 



1 

+ - 

P 



(2-6) 



Here the prime denotes derivative with respective to r]. 

For ^(77) = 1, the prepotential approach can generate exactly and quasi-exactly 
solvable systems associated with the classical orthogonal polynomials.® The pres- 
ence of ^ in the denominators of (p{x) and V{x) thus gives a rational extension, or 
deformation, of the traditional system. We therefore call ^(77) the deforming function. 

To make V exactly solvable, we demand that: (1) Wq is a regular function of x, 
(2) the deforming function ^{rj) has no zeros in the the ordinary (or physical) domain 
of ri(x), and (3) the function p{ri) does not appear in V. 

The requirement (3) can be easily met by setting the last term involving p{ri) in 
Eq. (|2-6p to a constant, say "— 



I.e., 



77V' + ( '^Wov + ii-2iQ'^^]p' + £p 



0. 



(2-7) 



If Wq, ^ can be determined, and Eq. (|2-7|) can be solved, then we would have con- 
structed an exactly solvable quantal system Hip = Sip defined hy H = —d?/dx'^ + 
V{x), with the wave function (|2-5p and the potential 



V{x) = Wl + 1^0 + 



2Wofj + i] 



(2.e 
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2.2. Determining Wo{x),^{i]) and p{rj) 

As mentioned before, we require that ^(r/) has no zeros in the ordinary domain 
of the variable 77, but ^(77) may have zeros in the other region in the complex ry-plane. 
Suppose ^{rj) satisfies the equation 



C2(^)^" + ci(7?)^' + 6'^ = 0, £ = real constant. 



(2- 



Here 02(7?) and ci(r/) are functions of ry to be determined. We want Eq. ()2-9p to be 
exactly solvable. This is most easily achieved by matching (|2-9|) with the (confluent) 
hypergeometric equation, and this we shall adopt in this paper. Thus C2(??) and ci(ry) 
are at most quadratic and linear in r/, respectively. 

If the factor exp(Wo(3^))/C(^) in Eq. ()2-5p is normalizable, then 75(77) = constant 
(in this case we shall take 75(17) = 1 for simplicity), which solves Eq. (|2-7p with f = 0, 
is admissible. This gives the ground state 



00 (x) 



(2-10) 



However, if exp(Wo{x)) / ^(rj) is non-normalizable, then (l)o{x) cannot be the ground 
state. In this case, the ground state, like all the excited states, must involve non- 
trivial 73(77) 7^ 1. Typically it is in such situation that the exceptional orthogonal 
polynomials arise. To determine states involving non-trivial p{rj), we proceed as 
follows. Since both ^ and ^' appear in Eq. (j2-7p . we shall take the ansatz that 75(77) 
be a linear combination of and 

p(7?) =e'(r7)F(77)+e(77)G(r/), (2-11) 

where -^(77) and G{r]) are some functions of rj. Then using Eq. ()2-9p we have 

£ 



ci 

C2 



F + F' + G 



-F + G' 



C2 



(2-12) 



We demand that Eq. (|2-7p be regular at the zeros of ^. This is achieved if p' oc ^, 
which requires that the coefficient of ^' in Eq. (j2-12p be zero, thus giving a relation 
that connects F and G, 



G='-^F 



F'. 



C2 



(2-13) 



Putting Eq. (12431) into 1^7}, we get 
.0 f £ 



■ 2 

^ drj 



C2 



£ 



-F + G' \ +£F 



-F + G'\+ (2Wof] + ii 




+ £G 



0. (2-14) 



Since and are independent for any polynomial ^, Eq. ()2T4p implies the coefficients 
of ^ and ^' are zero. Setting the terms in the two square-brackets to zero, and using 
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Eq. ()2-13p to eliminate G, we arrive at the following equations satisfied by F{ri) and 
ci(r/), respectively: 



r] F -\ ciF + 



and 



C2 



ci(^) = 

- El 



C2 



C2 



d / ci 



dr] \C2 



£F, 



(2-15) 



2dn 



2VFo?7 + V 
2Q{v) 



(2-16) 



where = Worj, and we have used the identity f/ = {drf / drj) /2 to arrive at the 
last line. 

Eq. ()2-13p suggests that we set 



F(r/) = C2{n)V{ri), 



(2-17) 



with some function V{ri) in order to avoid any possible singularity from C2- Eqs. (|2-15|) 
and ()2-13p then reduce to 



C2V" + (24 - ci) V + 



4 



V = 0, 



and 



Giv) 



4) V - C2V'. 



(2-18) 



(2-19) 



As mentioned before, in this paper we shall take €2(1]) and ci(r/) to be at most 
quadratic and linear in r/, respectively. This means the coefficients of the first and 
second terms in Eq. (j2-18p are also at most quadratic and linear in 77, respectively. So 
Eq. ()2T8p can be matched with the (confluent) hyper geometric equation, provided 
that the coefficient of the last term in ()2-18p is a constant. This then requires 



C2(??) 



(2-20) 



(Note: in general one has 02(7?) = ± constant x 7)^. But it is evident that the constant 
can be factored out of Eq. (|2T8p . together with a rescaling of £. Thus without loss 
of generality we set the constant to unity). From Eq. (j2-16p this leads to 



ciir]) 



1 d 
2dr] 



[if ] - 2Q(r?) 



(2-21) 



Now we summarize the procedure or algorithm for constructing an exactly solv- 
able quantal system, whose potential as well as its eigenfunctions and eigenvalues 
are all determined within the some approach: 

• choose 7]^ from a sinusoidal coordinate; this then fixes the form of C2; 
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by matching Eq. ()2-9p with the (confluent) hypergeometric equation, one deter- 
mines Q{r]), ci and ^(r/). Integrating Q{x) = WqT] then gives the prepotential 
Wq{x): 

Q{r}{x)) 



Wo{x) 



dx 

r][x) 



■2f V (2-22) 

• by matching Eq. ()2-18p with the (confluent) hypergeometric equation, one de- 
termines V, and thus F{r]), G{r]), p{r]) and £; 

• the exactly solvable system is defined by the wave function (|2-5p and the po- 
tential ()2-8p . which, by Eqs. (|2-9p and (|2-20p . can be recast in the form 



2VFo?? + ij]±ci 



±£. (2-23) 



2.3. Orthogonality of p{r]) 
Using the relations 



i dp . , (Pp .2 ,/,.., 

-7- = -, -i-=VP, -T^=riP +VP, 
arj fj ax ax^ 

one can recast Eq. (|2-7p into a differential equation in variable x, 



dx"^ 



p{r]{x)) + 2 [Wo 



S, I dx 



p{rj{x)) + £p{rj{x)) = 0. 



This can further be put in the Sturm-Liouville form 



d 

dx 



W2^p(r?(x)) 



where 



W{x) = exp 



+ SW^pivix)) = 0, 



dx\Wo-^ 



(2-24) 



(2-25) 



(2-26) 



According to the standard Sturm-Liouville theory, the functions ps (here we add a 
subscript to distinguish p corresponding to a particular eigenvalue £) are orthogonal. 



I.e., 



dx p£{r]{x))p£'{r]{x))W^{x) cx 5£^, 



in the x-space, or 
in the //-space. 



/ 



drj P£iv)P£'iv) 



W^{x{r])) 



(2-28) 
(2-29) 
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§3. L2 Laguerre case 

We now employ the above algorithm to generate the deformed radial oscillator 
as given in Ref . EJ-Hj) . 

Let us choose r/(x) = G [0, oo). Then t)^ = 4ry. For C2 and ci, we take 
the positive signs in Eqs. ()2-20p and ()2-2ip (the opposite situation is considered in 
Appendix B) . Thus C2(t/) = Ar] and ci = 2(1 — Q{rj)). 

3.1. Wq, i and £ 
Eq. (|2-9|) becomes 

^e" + ^(i-Q(^))e' + fe = o. (3-i) 

Comparing Eq. (|3-ip with the Laguerre equation 

7?4'(") + (a + l-r?)Lj")+^Lf) =0, £ = 0,1,2,..., (3-2) 
where L^^\r]) is the Laguerre polynomial, we have 

e(r?) = C,(r/;a) = Lf)(r/), £ = M, Q(ry) = 2 (^r/ - a - (3-3) 

For £,£{rj; a) not to have zeros in the ordinary domain [0, oo), we must have a < —i 
(see Appendix A). By Eq. ()2-22p . the form of Q{r]) gives 

x2 / 1 



Wo{x) = — -i^a + -jlnx. (3-4) 

We shall ignore the constant of integration as it can be absorbed into the normal- 
ization constant. 

3.2. pir]), 4>{fi) and £ 

The above results implies that exp(Wo) oc exp(x^/2)x~''"^2) (« < —t). The 
term exp(j;^/2) will make 0(x) non-normalizable if p{ri) = 1, or if V{ri) is a poly- 
nomial in r]. To remedy this, we try V = exp(—ri)U (rj) with some function U{rj). 
Eq. (|2-18p becomes 

rjU" + (-a + 1 - ry) [/' + (^^^ + U = 0. (3-5) 

Comparing this equation with the Laguerre equation ()3-2p (replacing i by another 
integer n = 0, 1, 2 . . .. Li the rest of this paper, the index n will always take on these 
values), one has 

U{r]) = 4-''\r]), £ = £n = ^n-a) -£ = 4.{n-a-e). (3-6) 
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From Eqs. ()2-17p and ()2-19p . one eventually obtains 

= Ae-'^Pe,n{r];a), (3-7) 

= r/4-)e; + («-n)4-°"'^e£- (3-8) 

Use has been made of Eqs. (|A-ip - (|A-3p in obtaining the last line in Eq. p-Sp . We 
note that Pi^niv'^ is a polynomial of degree i + n. It is just the L2 type exceptional 
Laguerre polynomial. We will show in the next subsection that it is equivalent to 
the form presented in Ref. [5]) (to be called HOS form for simplicity). By Eq. ()2-29p . 
one finds that Pi^niVj are orthogonal in the sense 

g-r?^-(Q+l) 

The exactly solvable potential is given by Eq. ()2-23p with Wo{x) and ^^(r/;a 
given by Eqs. ()3-4p and (|3-3p . respectively. The eigenvalues £„ are given in Eq. (j3-6p 
i.e. £n = 4(n — a — i) . Explicitly, the potential is 



dl] ^ P£,ni'n', a)Pf,^rri{m O) OC 5nm.- (3-9) 



+ 2(2^ -a). (3-10) 



It is easily shown that V{x) is equivalent to the potential for L2 Laguerre case in 
Ref. H]) ,[5]) . [TT]) with a = —g — (■ ~ \ (5 > 0) . Particularly, it is exactly equal to the 
form given in Eq. (2.21) of Ref. Ilip . The complete eigenfunctions are 

0^,„(x;a)oc P£,„(r/;a), a < -L (3-11) 

For ^ = 0, we have ^0 = 1 and = 0, and the system reduces to the radial 
oscillator. From Eq. ()3-8p one has P^^n Ln °^ and a < —1/2. 

3.3. Reducing Pi n{r]](^) to HOS form 

The polynomial Pi^ni^lj c^) is expressed as a bilinear combination of ^^(r/; a) and 
its derivative ^^(7/;q). The HOS form instead expresses the exceptional polynomial 
as a bilinear combination of ^£(??; a) and its shifted form, i.e., £,i{'r]; a — 1). 

To show the equivalence between P£^„(r/;a) and the HOS form, we make use of 
the identities Eqs. (|ATp and ()A-3p to express rj^iirj; a) in the first term of Pt^niVj 01) 
as 

rii'lima) = -r]Lf_X^\ri) 

= -aLg(r/)+£4"-')(77). (3-12) 



Then we have 



P,,„(r?;a) = (^-aLf\+ llf'^^) 4^") + (a^^) - r/L^")) L 



(a) 
I 



(a (4") - Lt\) + ILf-^^) 4""^ - r?4-")4"). (3-13) 



Prepotential approach to solvable rational potentials 



9 



Using Eq. ()A-2|1 we have L^^\r]) — L^^^{'q) = ^\fl)- Finally, we arrive at 

P,,„(r?;a) = (a + ^) 4-")(ry)e,(^; « " 1) " ^^^^""^(^KK^; «), (3-14) 

(3.15) 



Setting a = —g 



and ^e{r];g) = ^ '^'(v) ™to ()3-14p . we have 

1\ ^{ff+^+i), 



(3-16) 



This is, up to a multiplicative constant, the HOS form of the L2 Laguerre polynomial. 

The example in this section demonstrates that the prepotential approach de- 
scribed in Sect. 2 can indeed generate the exactly solvable quantal system which has 
the L2 Laguerre polynomials as the main part of its eigenfunctions. The prepoten- 
tial Wo{x), the potential V{x), the deforming function ^£{r];a), the eigenfunction 
(j)iAx;a) and eigenvalues f„ are all determined from first principle. 

In the next section and in the Appendix, we shall generate systems associated 
with the exceptional Jacobi and LI Laguerre polynomials. Our description for these 
cases will be concise, since the main steps are similar to those described in this 
section. 



§4. Exceptional Jacobi cases 



Let us choose rj{x) = cos(2x) G [—1, 1]. In this case it turns out, as can be easily 
confirmed, that both the upper and the lower signs in Eqs. ()2-20p and ()2-2ip for C2 
and ci give the same equations that determine ^ and V, i.e. Eqs. ()2-9p and (j2T8p . 
So for definiteness, we shall take the upper signs, which give C2(r/) = 4(1 — rj'^) and 
ci = -2(2r7 + Q(r?)). 

4.1. Wo, C and £ 

Equation determining ^ is 

(1 _ ^^)f^rj) + (-ry - ^) ^'(r?) + ^^v) = 0- (4-1) 

Comparing this with the differential equation satisfied by the Jacobi polynomial 
P^'^'^\r]), namely, 

(l-r?2)p;'(°''')(r/) + (/3-a-(a+/3+2)r/)P;("'^)(r?)+£(^+a+/3+l)pf •'')(?7) = 0, (4-2) 
we have 

av) = Uv, «, /5) = pt'^Hv), £ = m + « + /3 + 1), 

Q{ri) = 2 [a - /3 + (a + /3 + 1) 77] (4-3) 
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for some parameters a and /3. The form of Q{ri) gives, from Eq. ()2-22p . 

Wq{x) ~ ~ In sin X — ^/^ + In cos x. (4-4) 

The equation of V is 

(1 - 7?2)V" + [-/3 + a - (-/3 - a + 2) t?] V + [^-^ + « + /3 j V = 0. (4-5) 

From Eq. ()4-4p we have 

oc (1 - r/)-5("+l) (1 + . (4-6) 

The exponents in Eq. ()4-6p naturally divide the parameters a and j3 into four groups: 
(i) a > -1/2, /3 > -1/2, (ii) a > -1/2, ^ < -1/2, (iii) a < -1/2, /3 > -1/2 and 
(iv) a < —1/2, /3 < —1/2. Group (i) should be excluded, or ^i{r];a) will have zeros 
in the ordinary domain [—1, 1] (see Appendix A). So we shall study the other three 
cases. It turns out that these three cases correspond, respectively, to quantal systems 
related to the type Jl, J2 exceptional Jacobi polynomials, and a rationally extended 
Jacobi system obtained from the Darboux-Poschl- Teller system by deleting the lowest 
i excited states according to the Crum- Adler method discussed in Ref. [2T]) . We stress 
here that the actual admissible parameters a and /3 in each case are dictated by the 
final form of V, as will be shown below. 

4.2. Jl Jacobi case 

Consider the case with parameters a > —1/2, /? < —1/2. The deforming func- 
tion S^£{r];a, 13) is given in Eq. (j4-3p . 



Uv;a,f^) = Pt^\v)- (4-7) 

We demand that (,i{r];a, (3) has no zeros in the ordinary domain [—1,1]. From 
Eq. (|A-13p . one can easily check that this is the case if (3 < —i for a > —1/2. For 
this choice of the parameters the first term (1 — r/)" 2 ("+ 2) of Eq. (|4-6p will make the 
eigenfunction non-normalizable, if V is a polynomial. 

This prompted us to try V = (1 —rj)'^U{r]) where 7 is a real parameter and U{rj) 
a function of rj. From Eq. ()4-5p we find that U {rj) satisfies 

(1 - r/2)C/" + (-27 - /3 + a-(27-/?-a + 2)7?) U' 

+ (^^^ + a + /3 + 7(7 + /3-a-l) + 27(7 - 0)^3^^ U = 0. (4-8) 

If 7 = 0, Q, the coefficient of U in the last term of the above equation can be reduced 
to a constant, so that Eq. (j4-8p can be compared with the Jacobi diff^erential equation 
()4-2p . As 7 = does not solve our original problem with normalizability of the wave 
function, so we shall take 7 = a. This leads to 

(1 - r]^)U" + (-/3 - a - (-/3 + Q + 2)??) U' + { + /3(q + 1) W = 0. (4-9) 
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Comparing this equation with Eq. ()4-2p . we conclude that 

8=8n = ^ [n{n + a - /3 + 1) - + a + /3 + 1) - /3(a + 1)] . (4-10) 
Putting ah these results into F{r]) and G(r/) gives 

P,,„(7?;a,/5) ^ {(1 + r?)Pi"'-^)(r,)^; + [pPt^-^\r]) - (1 + r?)P^("'-«(r?)] C,} 

= (1 + r])Pt^-^Hr])ei - (n - /3))Pi"+^'-^-^)(^)e£. (4-11) 

We have made use of Eq. ()A-9p to arrive at the last line of ()4-lip . By Eq. ()2-29p . the 
orthogonality relations of P£^„(r/;a)'s are 

d^- k Pi,n{T,a,P)Pe,m{r];a,P) ^6nm. (4-12) 

1 Q 

The exactly solvable potential is given by Eq. (j2-23p with Wo{x) and ^i{r]]a,l3) 
given by Eqs. ()4-4p and ()4-3p . respectively. The eigenvalues f„ are given in Eq. ()4T0p . 
The complete eigenfunctions are 

<^,,„(x; a, /3) oc ^ '-^ ^—^ a, /3), (4-13) 

a > -1/2, /3 < -£. 

Using the identity (|A-9p one can show easily that 

Pe,n{v;a,(3) = {£ + (3)Pt'~^HvMmc^ + hf3-i) 

-(l + 7?)P^(-'-«(r/)eKr/;a,/3), (4-14) 

Up to a multiplicative constant, this is just the HOS form of the Jl Jacobi polynomial 
presented in Ref. [5]), with the substitution a = g + £ — 3/2 and P = —h — £ — ^. It 
is easy to show that V{x) and 8n are equivalent to those for Jl Jacobi case given in 
Ref. HI) ,[5]) ,[11]) with these values of a and /3. 

As ^ ^ 0, the system reduces to the trigonometric Darboux-Poschl- Teller po- 
tential, where a and f3 can now take the values a > —3/2, f3 < —1/2. 

4.3. J2 Jacobi case 

One can proceed in a similar manner to construct the exactly solvable systems 
with a < —1/2, /3 > —1/2. This turns out to lead to the system involving the J2 
Jacobi polynomials. 

We shall not bore the readers with similar details here. Instead, we point out that 
it is easier to obtain the system by symmetry consideration. One notes that under 
the parity transformation rj — >■ —r], together with interchange a -H- /3, Eqs. ()4-ip (with 
Q{ri) given by (j4-3p ) and (|4-5p are invariant in form. This is in complete accordance 
with the parity property of the Jacobi polynomials, namely 

p(-^l^)(-r,) = {-irpif^^-){r,). (4-15) 
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This implies that the J2 Jacobi system is simply the mirror image of the Jl Jacobi 
system, and thus it can be obtained from the Jl case by taking the above transfor- 
mations. 

4.4. Rationally extended Jacobi case 

Let a, /3 < —1/2. In this case, the factors in Eq. (|4-6p cause no problem with 
normalizability of the wave function even if V{ri) is a polynomial. For ^^(j?; a, (3) 
to be nodeless in the ordinary domain [—1, 1], we must choose a and j3 such that 
the conditions in ()A-13|) are satisfied. For example, if £ = 1, one can have a < 
-1, - 1 < /3 < -1/2, or /3 < -1, - 1 < a < -1/2. For ^ = 2, we have a, /3 < -2, 
or —2 < a, (3 < —1. For i odd, we must have a 7^ /?, or ^(r/) will have a zero at 
= in view of Eq. (j4-15p . 

Comparing Eqs. (|4-5|) and (|4-2|) . one obtains 

V(r?)=p(-°'-/^)(ry), 

£ = £^ = i [n(n -a-(3 + l)-iii + a + f3 + l)-a-(3]. (4-16) 

From F{ri) and G{r]) we get 

Piv) ^ ^^,n(r/;a,/3) = 4 {(1 - r?2)P^— /3)(r?)e; 
+ [(/3 - a - (/3 + a) r/) p(-°'-«(7?) - (1 - r/2)/^(-°'-«(^)] ^.j . (4-17) 
Again, by applying the identity (|A-7p and (|A-9p . one can reduce P£.„(r/; a, /?) to 
-P^,n(??;a,/3) 

= 4 {(£ + /3)(1 - r/)p(-"'-^)(r/)eKr/; a + 1, (3 - 1) 

+ (n-a)(l + r?)p(-"-^'-'^+i)(r/)e,(r/;a,/3)}. (4-18) 

One notes that Pi^niv'^ /5) is a polynomial of degree i + n + 1, and has n + 1 nodes. 
Thus the wave function with Pi^Q{i];a, (3) has one node, and does not correspond to 
the ground state. In fact, in this case the ground state wave function is given by 
Eq. ()2T0p with p{rj) = 1 and f = 0, since ^0(3;) is normalizable. To ensure that 
the energies of the excited states are positive, i.e., f „ > for n = 0,1,2,..., one 
must have, besides the constraints stated at the beginning of this subsection, the 
condition a + /3 < —i, which can be easily checked from the form of £ in Eq. ()4T6p . 

The functions Pf^„(r/; a, /?) (n = 0,1,2,...), together with p(r]) = 1, form a 
complete set and are orthogonal with respect to the weight function 

The complete eigenfunctions are given by 



(x; a, f3) oc 



it 



>(x; a, /3) cx ^ '-^ ^—^ a, /?). (4-20) 
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The exactly solvable potential is given by Eq. ()2-23p with Wo{x) and ^^(t/;q,/3) 
given by Eqs. (j4-4p and (j4-3p . respectively. Since the polynomials in the eigenfunc- 
tions start with degree zero, the polynomials P^^„(r/; a, /3) cannot be considered as 
exceptional. In fact, this system corresponds to the system discussed in Ref. [2T]) . 
which is obtained from the Darboux-Poschl- Teller system by deleting the lowest ^ 
excited states according to the Crum-Adler method. It belongs to the same class of 
rationally extended exactly solvable systems discussed in Ref. [22]) . [23]) . It is easy to 
show, using the identities in Appendix A. 2, that the polynomials in Eq. ()4T8p are 
proportional to those given by Eq. (A. 33) of Ref. 21) for the Jacobi case. 



We have demonstrated how all the quantal systems related to the exceptional 
Laguerre and Jacobi polynomials can be constructed in a direct and systematic way. 
In this approach one does not need to rely on the requirement of shape invariance 
and the Darboux-Crum transformation. Even the prepotential need not be assumed 
a priori. The prepotential, the deforming function, the potential, the eigenfunctions 
and eigenvalues are all derived within the same framework. It is worth to note that 
the main part of the eigenfunctions, which are the exceptional orthogonal polyno- 
mials, can be expressed as bilinear combination of the deformation function ^(77) 
and its derivative ^'{rf). However, they are equivalent of the forms given in Ref. [S]). 
We have also derived easily a rationally extended Jacobi model obtained from the 
Darboux-Poschl- Teller system by deleting the lowest ^ excited states according to 
the Crum-Adler method discussed in Rer. (2T|). 

We have not discussed the related hyperbolic Darboux-Poschl- Teller systems 
(which are of J2 type). They can be generated in the same way by choosing the 
appropriate sinusoidal coordinates. They can also be obtained from the trigonometric 
case by suitable analytic continuation. 



This work is supported in part by the National Science Council (NSC) of the 
Republic of China under Grant NSC NSC-99-2112-M-032-002-MY3. 



In this Appendix, we collect some useful identities satisfied by the Laguerre and 
Jacobi polynomials which are used in the main text. 

A.l. Laguerre Polynomials 

Some useful relations among Laguerre polynomials are: 



§5. Summary 
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Appendix A 

Useful identities 




(a+i) 



{ri), 



(A-l) 
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-1,-2, 



vLt\'\v) - (« + 1) Lt\'\ri) = M'^hv), (A-3) 

According to the Theorem 6.73 of Ref. for an arbitrary real number a ^ 
, the number of the positive zeros of L^"\ri) is £ if a > —1; it is 
^ + [a] + 1 if — £ < a < — 1; it is if a < —i. Here [a] denotes the integral part of a. 
Furthermore, = is a zero when and only when a = —1, —2, . . . , —i. 

A. 2. Jacobi polynomials 

Some useful relations among the Jacobi polynomial are: 

i + a + (3 + 1 



2(/3 + 1)P. 



~1 

2(£ + /3 + l)Pi"'''^(r,), 



(A-4) 



(A-5) 



he + a + (3 + l){v- l)Ptt''^^'\v), (A-6) 



(n + a)(l + r?)P,(!f +^)(r/) - /3(1 - r?)P;!+^''')(,?) = 2^P;"'^-^)(,?). 
Using Eqs. ()A-4p and ()A-6p to eliminate the P/"^"^''^"''"^ term gives 



(1 - r,) ^Pt''\v) = o.Pt'\^) + a)Pf -^'^^^^ 



dr] 



(A-7) 



(A-8) 



Combining Eqs. (jA-4p and (jA-Sp to eliminate the Pg^f'^'^\'i]) term gives 

(1 + ^Pt''^^'\v) = -(/? + l)P;"-''''+'^(^) + + /? + l)Pf'')(r?).(A.9) 



Setting a — >■ a + 1, /3 ^ /3 — 1, we get 



(1 + V) Y/t'\r]) = + /3)p("+^'^-^)(r?) - /3p("'«(ry). 



(A-10) 



According to the Theorem 6.72 of Ref. [Mp . for arbitrary real values of a and /3, 
the number of zeros of P^"'^\r]) in ( — 1, 1) is 



iVi(a,/3) 



2 [^] , if (-1)^ 
2 [f ] + 1, if (-1)^ 



£ 
£ 



£ + P 
£ 

£ + P 
£ 



>0; 



< 0. 



(A-ll) 



Here 



X = E 



(|2£ + a + /3 + l|-|a|- 1/31 + 1) 



(A-12) 
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where E{u) is the Klein's symbol defined by 

0, M < 0; 

E{u) = { [u] u > 0, u non-integral; 

u — 1 n = 1, 2, 3, . . . 

From this theorem, we conclude that the conditions for P^"'^\r]) to have no 
zeros in the ordinary domain ( — 1, 1) are 

|2£ + a + /3 + 1| - |a| - + 1 < 0, 

and(-l)'(^ + ")(^ + '')>0. (A.13) 

It is noted that 77 = +1(— 1) is a zero of P^°''^\r]) if and only if a(/3) = 
— 1, —2, ...,—£ with multiplicity \a\ (|/3|). 

Appendix B 

LI Laguerre case 



As with the L2 Laguerre case discussed in Sect. 3, let us take r]{x) = . But 
now the negative signs in Eqs. ()2-20p and (|2-2ip will be taken leading to C2(?y) = — 4r/ 
and ci = -2(1 - Q(r\)). 

B.l. Wo, i and £ 

Equation determining ^ is 

- vCiv) - ^ (1 - Qm eiv) + ^av) = o. (b-i) 

We shall take £ > 0, otherwise the problem reduces to the L2 Laguerre case discussed 
in Sect. 3. The first term of Eq. ()B-ip differs in sign from that of the the Laguerre 
equation (|3-2p . Suppose we make a parity change rj — )■ —rj in Eq. ()B-ip . then we will 
have 

VC"{-V) + ^ (1 - Q{-V)) a-v) + f = 0. (B-2) 

This equation has the form of Eq. (|3-2p , provided that 

a-ri)^U-iT,a) = Lf\n), £ = U, = 2 (77 - a - (B-3) 

for some parameter a. This means 

eKr/;a)=Lf)(-r/), (B-4) 

and 

Q(77) = -2 Tr? + a + (B-5) 



16 



C.-L. Ho 



For ^£(r/; a) not to have zeros in the ordinary domain [0, oo), we must have a > —1 
at the least (precise bound wih be determined later). The form of Qitj) then leads 
to 

Wo{x) = -^- (a + l)lnx. (B-6) 



2 V 2^ 
As before we ignore the constant of integration. 
B.2. p{r]), and 8 

Consider exp(Wo) oc exp(— a;2/2)x~''"~'"2) (a > —1 at the least). Contrary to 
the LI case, this time it is the term x~^"'*'2^ that could cause (t){x) non-normalizable 
(when a > —1/2) if V(r/) is a polynomial in 77. So we try V = r]^U{7]) where /? is a 
real parameter and f/(r/) a function of 77. Prom Eq. (|2-18p we get 

VU" + {2P-a + l-r,)U'+ + ^ _ /? _ ^7 = 0. (B-7) 

If /? = 0, a, the r/-dependent term in the last term of the above equation can be 
eliminated, and Eq. ()B-7p can be reduced to the Laguerre equation (|3-2p . As /3 = 
does not solve our original problem with normalizability of the wave function, we 
shall take /3 = a. This leads to 

?7(r/) = 4")(r/), 8 = £n = ^{n + a + l + l). (B-8) 

Putting all these result into F{rj) and G{rj) gives 

?>(??) =Pe,n{v) = -4r/"+^P^,„(??;a) 

p,,„(r?;a) = 4")e; + (L(r)-L;H)6 

= LWe; + 4"+^)6, (B.9) 

where use has been made of Eqs. (|A-ip and (|A-2p to get the last line. Pi^niw^^) 
is a polynomial of degree i + n. It will be shown below that it is just the LI 
type exceptional Laguerre polynomial. It is also easy to check that P£ „(r/;a)'s are 
orthogonal with respect to the weight function 

p-r; (a+l) 

-^J^- (B-10) 

The exactly solvable potential is given by Eq. ()2-23p with Wq{x) and ^£(r/;a) 
given by Eqs. (|B-6p and (|B-4p . respectively. The eigenvalues are £"„ = 4(n + a + £ + l) 
. It is easy to show that V{x) is equivalent to the potential for LI Laguerre case in 
Ref. H|l .[51) . [TT]) with a = g + I — 'i/2 {g > 0). Particularly, it is exactly equal to the 
form of potential in Eq. (2.20) of Ref. 11). The complete eigenfunctions are 

(j)i^n{x;a) oc P^,„(T?;a), a > --. (B-11) 

« 2 

As in the L2 case, this system reduces to the radial oscillator system in the limit 
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B.3. Reducing Pi^niv'j^) to HO S form 



Using Eqs. ()A-ip and (|A-2p . we have 



(B-12) 



Then it is easy to check that 



Pe,n{m «) = 4"^ ivMri; « + 1) - ivMm «) 



(B-13) 
(B.14) 



This is, up to a multiphcative constant, the HOS form of the LI Laguerre polynomial, 
with the substitution a = g + i — 3/2. 
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